Harmonic resist model for use in a lithographic apparatus and a device manufacturing method

ABSTRACT

A method for determining an image of a mask pattern in a resist coated on a substrate, the method including determining an aerial image of the mask pattern at substrate level; and convolving the aerial image with at least two orthogonal convolution kernels to determine a resist image that is representative of the mask pattern in the resist.

CROSS-REFERENCE TO RELATED APPLICATIONS

The present application is a Continuation of U.S. patent application Ser. No. 13/896,805, filed May 17, 2013 (now U.S. Pat. No. 8,682,059), which is a Divisional of U.S. patent application Ser. No. 12/625,079, filed Nov. 24, 2009, and issued as U.S. Pat. No. 8,447,095, which claims priority to U.S. Provisional Patent Application No. 61/117,396, filed Nov. 24, 2008, the contents of all of which are incorporated by reference in their entirety.

FIELD

The present invention generally relates to a method and computer product for performing simulations of the imaging results associated with a lithographic process.

BACKGROUND

A lithographic apparatus is a machine that applies a desired pattern onto a substrate, usually onto a target portion of the substrate. A lithographic apparatus can be used, for example, in the manufacture of integrated circuits (ICs). In such a case, a patterning device, which is alternatively referred to as a mask or a reticle, may be used to generate a circuit pattern to be formed on an individual layer of the IC. This pattern can be transferred onto a target portion (e.g. including part of, one, or several dies) on a substrate (e.g. a silicon wafer). Transfer of the pattern is typically via imaging onto a layer of radiation-sensitive material (resist) provided on the substrate. In general, a single substrate will contain a network of adjacent target portions that are successively patterned. Conventional lithographic apparatus include so-called steppers, in which each target portion is irradiated by exposing an entire pattern onto the target portion at once, and so-called scanners, in which each target portion is irradiated by scanning the pattern through a radiation beam in a given direction (the “scanning”-direction) while synchronously scanning the substrate parallel or anti-parallel to this direction. It is also possible to transfer the pattern from the patterning device to the substrate by imprinting the pattern onto the substrate.

In order to develop a practical lithographic resist process, it is desirable to first use a lithography simulation tool for process optimization. Subsequently, actual substrate printing verification may be performed. The challenge is how to ensure the prediction of a resist critical dimension (CD), a contour of the resist pattern and/or whether some contours of resist patterns are merged or cut with sufficient accuracy when using the simulation tool. Resist model calibration is one of the critical key factors in the simulation process, and the robustness of the simulation model is another important factor.

A typical lithography simulation involves three basic steps. First, an aerial image for the feature in question is calculated. The aerial image calculation is based on the optical settings of a lithographic apparatus, which include, for example, numerical aperture, exposure wavelength and characteristics of lenses. Second, a post exposure bake (PEB) step is performed. This step provides two functions: (1) allowing chemical amplification for photo speed to take place due to the heat, and (2) minimizing resist CD swing caused by standing wave effects. The third step involves developing a resist pattern based on the diffused aerial image.

In general, there are two conventional approaches to resist modeling: either through a slow but more physical modeling of the process or through a faster but empirical approach. However, none of these approaches to resist modeling provides satisfactory results. Indeed, these approaches may not fully represent the chemical processes and/or are generally time-consuming.

SUMMARY

It is desirable to provide a method that accurately and quickly simulates the image of a mask pattern in a photoresist.

In an aspect of the invention, there is provided a method for determining an image of a mask pattern in a resist coated on a substrate, the method including determining an aerial image of the mask pattern at substrate level; and convolving the aerial image with at least two orthogonal convolution kernels to determine a resist image that is representative of the mask pattern in the resist.

In another aspect of the invention, there is provided a method for determining an image of a mask pattern in a resist coated on a substrate, the method including determining an aerial image of the mask pattern at substrate level; and transforming the aerial image into a resist image that is representative of the image of the mask pattern in the resist, the aerial image being transformed into the resist image using a model that includes at least two convolution kernels, the transforming including convolving the aerial image with the at least two convolution kernels, wherein the at least two convolution kernels are selected such that the transformation of the aerial image into the resist image has the properties of rotational and mirror symmetry conservation.

In another aspect of the invention, there is provided a computer product having machine executable instructions, the instructions being executable by a machine to perform a method for determining an image of a mask pattern in a resist coated on a substrate, the method including determining an aerial image of the mask pattern at substrate level; and convolving the aerial image with at least two orthogonal convolution kernels to determine a resist image that is representative of the image of the mask pattern in the resist.

In another embodiment of the invention, there is provided a computer product having machine executable instructions, the instructions being executable by a machine to perform a method for determining an image of a mask pattern in a resist coated on a substrate, the method including determining an aerial image of the mask pattern at substrate level; and transforming the aerial image into a resist image that is representative of the image of the mask pattern in the resist, the aerial image being transformed into the resist image using a model that includes at least two convolution kernels, the transforming including convolving the aerial image with the at least two convolution kernels, wherein the at least two convolution kernels are selected such that the transformation of the aerial image into the resist image has the properties of rotational and mirror symmetry conservation.

These and other aspects, features, and characteristics of the present invention, as well as the methods of operation and functions of the related elements of structure will become more apparent upon consideration of the following description and the appended claims with reference to the accompanying drawings, all of which form a part of this specification, wherein like reference numerals designate corresponding parts in the various figures. It is to be expressly understood, however, that the drawings are for the purpose of illustration and description only and are not intended as a definition of the limits of the invention. As used in the specification and in the claims, the singular form of “a”, “an”, and “the” include plural referents unless the context clearly dictates otherwise.

BRIEF DESCRIPTION OF THE DRAWINGS

Embodiments of the invention will now be described, by way of example only, with reference to the accompanying schematic drawings in which corresponding reference symbols indicate corresponding parts, and in which:

FIG. 1 shows a lithographic apparatus in accordance with an embodiment of the invention;

FIG. 2 is an exemplary flowchart illustrating a process of obtaining a resist image in accordance with an embodiment of the invention;

FIG. 3 illustrates the rotational property of the process of FIG. 2;

FIG. 4 illustrates the mirror symmetry property of the process of FIG. 2; and

FIG. 5 is an exemplary flowchart illustrating a process of obtaining a resist image in accordance with an embodiment of the invention.

DETAILED DESCRIPTION

FIG. 1 schematically depicts a lithographic apparatus according to one embodiment of the invention. The apparatus includes an illumination system (illuminator) IL configured to condition a radiation beam B (e.g. UV radiation or any other suitable radiation), a patterning device support or mask support structure (e.g. a mask table) MT constructed to support a patterning device (e.g. a mask) MA and connected to a first positioning device PM configured to accurately position the patterning device in accordance with certain parameters. The apparatus also includes a substrate table (e.g. a wafer table) WT or “substrate support” constructed to hold a substrate (e.g. a resist-coated wafer) W and connected to a second positioning device PW configured to accurately position the substrate in accordance with certain parameters. The apparatus further includes a projection system (e.g. a refractive projection lens system) PS configured to project a pattern imparted to the radiation beam B by patterning device MA onto a target portion C (e.g. including one or more dies) of the substrate W.

The illumination system may include various types of optical components, such as refractive, reflective, magnetic, electromagnetic, electrostatic or other types of optical components, or any combination thereof, to direct, shape, or control radiation.

The patterning device support holds the patterning device in a manner that depends on the orientation of the patterning device, the design of the lithographic apparatus, and other conditions, such as for example whether or not the patterning device is held in a vacuum environment. The patterning device support can use mechanical, vacuum, electrostatic or other clamping techniques to hold the patterning device. The patterning device support may be a frame or a table, for example, which may be fixed or movable as required. The patterning device support may ensure that the patterning device is at a desired position, for example with respect to the projection system. Any use of the terms “reticle” or “mask” herein may be considered synonymous with the more general term “patterning device.”

The term “patterning device” used herein should be broadly interpreted as referring to any device that can be used to impart a radiation beam with a pattern in its cross-section so as to create a pattern in a target portion of the substrate. It should be noted that the pattern imparted to the radiation beam may not correspond exactly to the desired pattern in the target portion of the substrate, for example if the pattern includes phase-shifting features or so called assist features. Generally, the pattern imparted to the radiation beam will correspond to a particular functional layer in a device being created in the target portion, such as an integrated circuit.

The patterning device may be transmissive or reflective. Examples of patterning devices include masks, programmable mirror arrays, and programmable LCD panels. Masks are well known in lithography, and include mask types such as binary, alternating phase-shift, and attenuated phase-shift, as well as various hybrid mask types. An example of a programmable mirror array employs a matrix arrangement of small mirrors, each of which can be individually tilted so as to reflect an incoming radiation beam in different directions. The tilted mirrors impart a pattern in a radiation beam which is reflected by the mirror matrix.

The term “projection system” used herein should be broadly interpreted as encompassing any type of projection system, including refractive, reflective, catadioptric, magnetic, electromagnetic and electrostatic optical systems, or any combination thereof, as appropriate for the exposure radiation being used, or for other factors such as the use of an immersion liquid or the use of a vacuum. Any use of the term “projection lens” herein may be considered as synonymous with the more general term “projection system”.

As depicted here, the apparatus is of a transmissive type (e.g. employing a transmissive mask). Alternatively, the apparatus may be of a reflective type (e.g. employing a programmable mirror array of a type as referred to above, or employing a reflective mask).

The lithographic apparatus may be of a type having two (dual stage) or more substrate tables or “substrate supports” (and/or two or more mask tables or “mask supports”). In such “multiple stage” machines the additional tables or supports may be used in parallel, or preparatory steps may be carried out on one or more tables or supports while one or more other tables or supports are being used for exposure.

The lithographic apparatus may also be of a type wherein at least a portion of the substrate may be covered by a liquid having a relatively high refractive index, e.g. water, so as to fill a space between the projection system and the substrate. An immersion liquid may also be applied to other spaces in the lithographic apparatus, for example, between the patterning device (e.g. mask) and the projection system Immersion techniques can be used to increase the numerical aperture of projection systems. The term “immersion” as used herein does not mean that a structure, such as a substrate, must be submerged in liquid, but rather only means that a liquid is located between the projection system and the substrate during exposure.

Referring to FIG. 1, the illuminator IL receives a radiation beam from a radiation source SO. The source and the lithographic apparatus may be separate entities, for example when the source is an excimer laser. In such cases, the source is not considered to form part of the lithographic apparatus and the radiation beam is passed from the source SO to the illuminator IL with the aid of a beam delivery system BD including, for example, suitable directing mirrors and/or a beam expander. In other cases the source may be an integral part of the lithographic apparatus, for example when the source is a mercury lamp. The source SO and the illuminator IL, together with the beam delivery system BD if required, may be referred to as a radiation system.

The illuminator IL may include an adjuster AD configured to adjust the angular intensity distribution of the radiation beam. Generally, at least the outer and/or inner radial extent (commonly referred to as σ-outer and σ-inner, respectively) of the intensity distribution in a pupil plane of the illuminator can be adjusted. In addition, the illuminator IL may include various other components, such as an integrator IN and a condenser CO. The illuminator may be used to condition the radiation beam to have a desired uniformity and intensity distribution in its cross-section.

The radiation beam B is incident on the patterning device (e.g., mask) MA, which is held on the patterning device support (e.g., mask table) MT, and is patterned by the patterning device. Having traversed the patterning device (e.g. mask) MA, the radiation beam B passes through the projection system PS, which focuses the beam onto a target portion C of the substrate W. With the aid of the second positioning device PW and position sensor IF (e.g. an interferometric device, linear encoder or capacitive sensor), the substrate table WT can be moved accurately, e.g. so as to position different target portions C in the path of the radiation beam B. Similarly, the first positioning device PM and another position sensor (which is not explicitly depicted in FIG. 1) can be used to accurately position the patterning device (e.g. mask) MA with respect to the path of the radiation beam B, e.g. after mechanical retrieval from a mask library, or during a scan. In general, movement of the patterning device support (e.g. mask table) MT may be realized with the aid of a long-stroke module (coarse positioning) and a short-stroke module (fine positioning), which form part of the first positioning device PM. Similarly, movement of the substrate table WT or “substrate support” may be realized using a long-stroke module and a short-stroke module, which form part of the second positioner PW. In the case of a stepper (as opposed to a scanner) the patterning device support (e.g. mask table) MT may be connected to a short-stroke actuator only, or may be fixed. Patterning device (e.g. mask) MA and substrate W may be aligned using patterning device alignment marks M1, M2 and substrate alignment marks P1, P2. Although the substrate alignment marks as illustrated occupy dedicated target portions, they may be located in spaces between target portions (these are known as scribe-lane alignment marks). Similarly, in situations in which more than one die is provided on the patterning device (e.g. mask) MA, the patterning device alignment marks may be located between the dies.

The depicted apparatus could be used in at least one of the following modes:

1. In step mode, the patterning device support (e.g. mask table) MT or “mask support” and the substrate table WT or “substrate support” are kept essentially stationary, while an entire pattern imparted to the radiation beam is projected onto a target portion C at one time (i.e. a single static exposure). The substrate table WT or “substrate support” is then shifted in the X and/or Y direction so that a different target portion C can be exposed. In step mode, the maximum size of the exposure field limits the size of the target portion C imaged in a single static exposure.

2. In scan mode, the patterning device support (e.g. mask table) MT or “mask support” and the substrate table WT or “substrate support” are scanned synchronously while a pattern imparted to the radiation beam is projected onto a target portion C (i.e. a single dynamic exposure). The velocity and direction of the substrate table WT or “substrate support” relative to the patterning device support (e.g. mask table) MT or “mask support” may be determined by the (de-)magnification and image reversal characteristics of the projection system PS. In scan mode, the maximum size of the exposure field limits the width (in the non-scanning direction) of the target portion in a single dynamic exposure, whereas the length of the scanning motion determines the height (in the scanning direction) of the target portion.

3. In another mode, the patterning device support (e.g. mask table) MT or “mask support” is kept essentially stationary holding a programmable patterning device, and the substrate table WT or “substrate support” is moved or scanned while a pattern imparted to the radiation beam is projected onto a target portion C. In this mode, generally a pulsed radiation source is employed and the programmable patterning device is updated as required after each movement of the substrate table WT or “substrate support” or in between successive radiation pulses during a scan. This mode of operation can be readily applied to maskless lithography that utilizes programmable patterning devices, such as a programmable mirror array of a type as referred to above.

Combinations and/or variations on the above described modes of use or entirely different modes of use may also be employed.

Conventional approaches that determine the resist image in a lithographic simulation involve convolving the aerial image (or a pixel-by-pixel transform of the aerial image) with various weighted kernels. These weighted kernels are usually Gaussian functions. These approaches are generally undesirable for several reasons. First, the chemical processes of conventional models are generally not computed accurately. With that said, the accuracy of the models cannot be increased without significantly increasing the computation time, which is undesirable. As a result, an empirical model is generally used. Second, the fitting of conventional empirical models is generally difficult. This results in potentially unstable fitting procedures, characterized by cancellations of model terms. Current models can have two Gaussian kernels with similar standard deviations. Two such kernels will be highly correlated and numerically unstable. Third, the kernels are selected arbitrarily and do not have any physical meaning. Fourth, current terms of simulation models have many non-linear parameters that need to be fitted. This requires many interactions between the user and the computer program. The user may have to run several fitting processes, changing the search range of the non-linear terms and verifying that the results are accurate and physically realistic. This approach is tedious and time-consuming for the user. This approach further requires more training and it is more an art than a science. The results of the model calibration may be highly dependent on who is doing it, which is undesirable in a production environment.

In one aspect of the invention, there is provided a model that accurately and quickly simulates the image of a mask pattern in a photoresist. The model, according to one embodiment of the invention, does not refer to any chemical processes. Rather, it is based mainly on mathematical and symmetry considerations of the underlying physics. The model involves convolving an aerial image with convolution kernels that form a complete orthogonal basis. It will be appreciated that a basis is complete when an infinite number of kernels are present. In one embodiment, some of the kernels may be used to calculate the resist image. The use of convolution orthogonal kernels to simulate a resist image provides significant benefits and unexpected results in terms of computation speed, accuracy and user's interaction over conventional simulation models.

Referring now to FIG. 2, this figure shows an exemplary flowchart illustrating a process 200 of obtaining a resist image of a mask pattern in accordance with an embodiment of the invention. The process 200 begins at procedure 205 where an aerial image of a mask pattern is determined at substrate level. The determined aerial image consists of a bitmap image of the pattern. The bitmap image is representative of the light intensity distribution of the image of the mask pattern projected by the lithographic apparatus. In practice, the aerial image of the mask pattern can be calculated using known lithography software simulation technologies/programs, such as Brion's Tachyon™, or Lithocruiser™, or KLA-Tencor's Prolith™. The contours of the mask pattern to be imaged, which may include optical proximity corrections (OPC), is provided in a standardized digital file format such as GDS, GDSII or OASIS and is inputted into these lithography software simulation programs to calculate the aerial image.

Determination of the aerial image of the mask pattern by these software products involves mathematically modeling the imaging of the pattern onto the substrate (e.g. from a mask), taking into account the illumination source, the projection optics and the pattern details. Efficient simulation methods that are fast enough for full-chip applications approximate the realistic 3-dimensional intensity distribution on the substrate by a 2-dimensional aerial image. An efficient implementation of a lithography model is possible using the following formalism, where the image (here in scalar form, which may be extended to include polarization vector effects) is expressed as a Fourier sum over signal amplitudes in the pupil plane. According to the standard Hopkins theory, the aerial image may be defined by:

$\begin{matrix} \begin{matrix} {{A(x)} = {\sum_{k}{{S(k)}{{\sum_{k^{\prime}}{{M\left( {k^{\prime} - k} \right)}{P\left( k^{\prime} \right)}{\exp\left( {{- j}\; k^{\prime}x} \right)}}}}^{2}}}} \\ {= {\sum_{k}{{S(k)}\left\{ {\sum_{k^{\prime}}{\sum_{k^{''}}{{M\left( {k^{\prime} - k} \right)}{P\left( k^{\prime} \right)}{M^{*}\left( {k^{''} - k} \right)}{P^{*}\left( k^{''} \right)}{\exp\left( {{- {j\left( {k^{\prime} - k^{''}} \right)}}x} \right)}}}} \right\}}}} \\ {= {\sum_{k^{\prime}}{\sum_{k^{''}}\left\lbrack {\sum_{k}{{S(k)}{P\left( {k + k^{\prime}} \right)}{P^{*}\left( {k + k^{''}} \right)}}} \right\rbrack}}} \\ {{M\left( k^{\prime} \right)}{M^{*}\left( k^{''} \right)}{\exp\left( {{- {j\left( {k^{\prime} - k^{''}} \right)}}x} \right)}} \\ {= {\sum_{k^{\prime}}{\sum_{k^{''}}{{TCC}_{k^{\prime},k^{''}}{M\left( k^{\prime} \right)}{M^{*}\left( k^{''} \right)}}}}} \\ {\exp\left( {{- {j\left( {k^{\prime} - k^{''}} \right)}}x} \right)} \end{matrix} & \left( {{Eq}.\mspace{14mu} 1} \right) \end{matrix}$ where, A (x) is the aerial image intensity at point x within the image plane (for notational simplicity, a two-dimensional coordinate represented by a single variable is utilized), k represents a point on the source plane, S (k) is the source intensity from point k, k′ and k″ are points on the pupil plane, M is the Fourier transform of the mask image, P is the pupil function, and TCC_(k′,k″)=Σ_(k) S(k) P(k+k′) P*(k+k″). An important aspect of the foregoing derivation is the change of summation order (moving the sum over k inside) and indices (replacing k′ with k+k′ and replacing k″ with k+k″), which results in the separation of the Transmission Cross Coefficients (TCCs), defined by the term inside the square brackets in the third line of the equation. These coefficients are independent of the mask pattern and therefore can be pre-computed using knowledge of the optical elements or configuration only (e.g., NA and σ or the detailed illuminator profile). It is further noted that although in the given example (Eq. 1) is derived from a scalar imaging model, this formalism can also be extended to a vector imaging model, where the polarized light components are summed separately.

Furthermore, the approximate aerial image can be calculated by using only a limited number of dominant TCC terms, which can be determined by diagonalizing the TCC matrix and retaining the terms corresponding to its largest eigenvalues, i.e.,

$\begin{matrix} {{TCC}_{k^{\prime},k^{''}} = {\sum\limits_{i = 1}^{N}{\lambda_{i}{\phi_{i}\left( k^{\prime} \right)}{\phi_{i}^{*}\left( k^{''} \right)}}}} & \left( {{Eq}.\mspace{14mu} 2} \right) \end{matrix}$ where λ_(i) (i=1, . . . , N) denotes the N largest eigenvalues and φ_(i)(.) denotes the corresponding eigenvector of the TCC matrix. It is noted that (Eq. 2) is exact when all terms are retained in the Eigen series expansion, i.e., when N is equal to the rank of the TCC matrix. However, in actual applications, it is typical to truncate the series by selecting a smaller N to increase the speed of the computation process.

$\begin{matrix} {{A(x)} = {\sum\limits_{k^{\prime}}{\sum\limits_{k^{''}}{{TCC}_{k^{\prime},k^{''}}{M\left( k^{\prime} \right)}{M^{*}\left( k^{''} \right)}{\exp\left( {{- {j\left( {k^{\prime} - k^{''}} \right)}}x} \right)}}}}} \\ {= {\sum\limits_{k^{\prime}}{\sum\limits_{k^{''}}{\sum\limits_{i = 1}^{N}{\lambda_{i}{\phi_{i}\left( k^{\prime} \right)}{\phi_{i}^{*}\left( k^{''} \right)}{M\left( k^{\prime} \right)}{M^{*}\left( k^{''} \right)}{\exp\left( {{- {j\left( {k^{\prime} - k^{''}} \right)}}x} \right)}}}}}} \\ {= {\sum\limits_{i = 1}^{N}{\lambda_{i}{\sum\limits_{k^{\prime}}{{\phi_{i}\left( k^{\prime} \right)}{M\left( k^{\prime} \right)}{\exp\left( {{- j}\; k^{\prime}x} \right)}}}}}} \\ {\sum\limits_{k^{''}}{{\phi_{i}^{*}\left( k^{''} \right)}{M^{*}\left( k^{''} \right)}{\exp\left( {j\; k^{''}x} \right)}}} \\ {= {\sum\limits_{i = 1}^{N}{\lambda_{i}{{\Phi_{i}(x)}}^{2}}}} \end{matrix}$

Thus, (Eq. 1) can be rewritten as:

(Eq. 3)

where

${\Phi_{i}(x)} = {\sum\limits_{k^{''}}{{\phi_{i}\left( k^{''} \right)}{M\left( k^{''} \right)}{\exp\left( {{- j}\; k^{''}x} \right)}}}$ and |.| denotes the magnitude of a complex number.

Using a sufficiently large number of TCC terms and a suitable model calibration methodology allows for an accurate description of the optical projection process and provides an accurate representation of the aerial image. Additional information regarding the calculation of an aerial image can be gleaned from U.S. Pat. No. 6,871,337 issued on Mar. 22, 2005 titled “Illumination Optimization For Specific Mask Patterns”, the contents of which is incorporated herein in its entirety by reference.

Referring back to FIG. 2, after determining the aerial image A of the mask pattern, the method proceeds to procedure 210 where the aerial image A is convolved with a number of orthogonal convolution kernels (e.g. in an embodiment at least two orthogonal convolution kernels) to determine a resist image R that is representative of the pattern in the resist. As discussed above, the aerial image A of the mask pattern is a bitmap image that represents the light intensity distribution of the image of the mask pattern projected by the lithographic apparatus. Convolution of each point of the bitmap image with a convolution kernel operates to smooth the energy of that point within the resist and provides a second bitmap image that is representative of the image of the mask pattern in the resist. In essence, the convolution of the aerial image of the mask pattern with a convolution kernel results in averaging the intensity of energy in the resist—averaging the energy at a point and energies nearby.

This resist image can then be used to perform optical proximity correction, contour computation, defect detection, or verification of manufacturability.

In an embodiment of the invention, the convolution kernels that are used to convolve the aerial image are orthogonal kernels that form an orthogonal basis. In one implementation, these orthogonal convolution kernels are solutions ψ of a two dimensional quantum harmonic oscillator ĤΨ=EΨ, with

$\hat{H} = {{{- \frac{\hslash^{2}}{2m}}\Delta_{2}} + {\frac{1}{2}{mw}^{2}{r^{2} \cdot \left( {\Delta_{2} = {{\frac{1}{r}\frac{\partial}{\partial r}} + \frac{\partial^{2}}{\partial r^{2}} + {\frac{1}{r^{2}}\frac{\partial^{2}}{\partial\theta^{2}}}}} \right)}}}$

The solutions ψ can be defined in polar coordinates as Ψ(r,θ)=R(r)Θ(θ).

$\begin{matrix} {\frac{\Theta^{''}}{\Theta} = {{{{- \frac{2m}{\hslash^{2}}}\left( {E - {\frac{1}{2}{mw}^{2}r^{2}}} \right)r^{2}} - {r\;\frac{R^{\prime}}{R}} - {r^{2}\;\frac{R^{''}}{R}}} = {{constant} = {{{{- {n^{2}.\left( {{{Let}\mspace{14mu}\Theta} = {\mathbb{e}}^{i\; n\;\Theta}} \right)}}\mspace{14mu}{or}\mspace{14mu} R^{''}} + {\frac{1}{r}R^{\prime}} + {\frac{2m}{h^{2}}\left( {E - {\frac{1}{2}{mw}^{2}r^{2}} - \frac{n^{2}}{r^{2}}} \right)R}} = 0}}}} & \left( {{Eq}.\mspace{14mu} 4} \right) \end{matrix}$

Assuming that

$\alpha \equiv \frac{mw}{h}$ and R takes the form of

${R = {{\exp\left( {{- \frac{\alpha}{2}}r^{2}} \right)}{\mathcal{R}(r)}}},$ equation 4 becomes

${\mathcal{R}^{''} + {\frac{1}{r}\left( {1 - {2\alpha\; r^{2}}} \right)\mathcal{R}^{\prime}} + {\left( {\frac{2{mE}}{h^{2}} - {2\alpha} - \frac{n^{2}}{r^{2}}} \right)\mathcal{R}}} = 0$

It will be appreciated that the radius r can be scaled by the factor r₀, with r₀=√{square root over (2/α)}, in which case

$R = {{\exp\left( {- \frac{r^{2}}{r_{0}^{2}}} \right)}{{\mathcal{R}\left( \frac{r}{r_{0}} \right)}.}}$

Assuming

is a function of αr²:x≡αr²,

≡y(x), the function above becomes

$\begin{matrix} {{{xy}^{''} + {\left( {1 - x} \right)y^{\prime}} + {\left( {\frac{E}{2w\;\hslash} - \frac{1}{2} - \frac{n^{2}}{4x}} \right)y}} = 0} & \left( {{Eq}.\mspace{14mu} 5} \right) \end{matrix}$ When n=0, this equation is the Laguerre's equation, and has the Laguerre polynomials as its solution.

Equation 5 can be rewritten in the following manner:

$\begin{matrix} {{{xy}^{''} + {\left( {1 - x} \right)y^{\prime}} + {\left( {E_{k} - \frac{n^{2}}{x}} \right)y}} = 0} & \left( {{Eq}.\mspace{14mu} 6} \right) \end{matrix}$

When n=0, the solutions of equation 6 are Laguerre polynomials L_(k)(x). When n≠0, the solutions are based on the Whittaker functions:

$\frac{{\mathbb{e}}^{x/2}W_{{k + {1/2}},{n{(x)}}}}{\sqrt{(x)}}$

The Whittaker function is given by:

$\begin{matrix} {{W_{{k + {1/2}},n}(x)} = {x^{{1/2} + n}{{\mathbb{e}}^{{- x}/2}\left\lbrack {1 + {\frac{n - k}{{1!}\left( {{2n} + 1} \right)}x} + {\frac{\left( {n - k} \right)\left( {n - k + 1} \right)}{{2!}\left( {{2n} + 1} \right)\left( {{2n} + 2} \right)}x^{2}} + \ldots}\mspace{14mu} \right\rbrack}}} & \left( {{Eq}.\mspace{14mu} 7} \right) \end{matrix}$

This simplifies to the solution S_(k,n)(x) as follows:

$\begin{matrix} {{S_{k,n}(x)} = {x^{n}\left\lbrack {1 + {\frac{n - k}{{1!}\left( {{2n} + 1} \right)}x} + {\frac{\left( {n - k} \right)\left( {n - k + 1} \right)}{{2!}\left( {{2n} + 1} \right)\left( {{2n} + 2} \right)}x^{2}} + \ldots}\mspace{14mu} \right\rbrack}} & \left( {{Eq}.\mspace{14mu} 8} \right) \end{matrix}$ When n=0, the Laguerre polynomials are recovered. The polynomials are finite.

The solution S_(k,n)(x) are orthogonal given an exponential weight function:

$\begin{matrix} {{{\int_{0}^{\infty}{{\mathbb{e}}^{- x}{S_{k,n}(x)}{S_{k^{\prime},n^{\prime}}(x)}{\mathbb{d}x}}} = {{\delta\left\lbrack {k - k^{\prime}} \right\rbrack}{\delta\left\lbrack {n - n^{\prime}} \right\rbrack}}},} & \left( {{Eq}.\mspace{14mu} 9} \right) \end{matrix}$ where k and n are integers and δ is a discrete Dirac indicator, δ[0]=1 and δ[i]=0 for i≠0.

This property is used to create two dimensional kernels K_(k,n)(r,θ) in accordance with an embodiment of the invention. These kernels are defined as follows: K _(k,n)(r,θ)^(±) =e ^(−r) ² ^(/2) S _(k,n)(r ²)e ^(±inθ)

The orthogonality of the kernels K_(k,n)(r,θ) can be verified in the following manner:

$\begin{matrix} {{\int_{- \infty}^{\infty}{\int_{- \infty}^{\infty}{{K_{k,n}\left( {x,y} \right)}{K_{k^{\prime},n^{\prime}}\left( {x,y} \right)}{\mathbb{d}x}{\mathbb{d}y}}}} = {\int_{0}^{\infty}{\int_{0}^{2\Pi}{{K_{k,n}\left( {r,\theta} \right)}{K_{k^{\prime},n^{\prime}}\left( {r,\theta} \right)}r{\mathbb{d}r}{\mathbb{d}\theta}}}}} \\ {= {\int_{0}^{\infty}{{\mathbb{e}}^{{- r^{2}}/2}{S_{k,n}\left( r^{2} \right)}{\mathbb{e}}^{{- r^{2}}/2}{S_{k^{\prime},n^{\prime}}\left( r^{2} \right)}r{{\mathbb{d}r} \cdot}}}} \\ {\int_{0}^{2\Pi}{{\mathbb{e}}^{{\pm {in}}\;\theta}{\mathbb{e}}^{{\pm {in}^{\prime}}\theta}{\mathbb{d}\theta}}} \\ {= {\frac{1}{2}{\int_{0}^{\infty}{{\mathbb{e}}^{- u}{S_{k,n}(u)}{S_{k^{\prime},n^{\prime}}(u)}{{\mathbb{d}u} \cdot}}}}} \\ {\int_{0}^{2\Pi}{{\mathbb{e}}^{{\pm {in}}\;\theta}{\mathbb{e}}^{{\pm {in}^{\prime}}\theta}{\mathbb{d}\theta}}} \\ {= {{\delta\left\lbrack {k - k^{\prime}} \right\rbrack}{\delta\left\lbrack {n - n^{\prime}} \right\rbrack}}} \end{matrix}$ When n=0, the solutions (based on Laguerre polynomials) form a complete basis of all the rotationally independent functions.

The aerial image A is convolved with one or more convolution kernels K_(k,n)(r,θ) to determine the resist image R, i.e. the image of the mask pattern in the resist. In an embodiment, the resist image R is calculated in the following manner:

$\begin{matrix} {R = {{\sum\limits_{i = 0}^{\infty}{c_{i}\left( {A*K_{i}} \right)}} + {\sum\limits_{i = 0}^{\infty}{{d_{i}\left( {A*K_{1i}} \right)} \cdot \left( {A*K_{2i}} \right)}}}} & \left( {{Eq}.\mspace{14mu} 10} \right) \end{matrix}$ wherein R is a bitmap resist image, A is a bitmap aerial image of the mask pattern at substrate level, K_(i), K_(1i), K_(2i) are orthogonal convolution kernels and c_(i) and d_(i) are fitting coefficients, where 0≦i≦∞. The first sum

$\sum\limits_{i = 0}^{\infty}{c_{i}\left( {A*K_{i}} \right)}$ is referred to as linear, while the second sum

$\sum\limits_{i = 0}^{\infty}{{d_{i}\left( {A*K_{1i}} \right)} \cdot \left( {A*K_{2i}} \right)}$ is referred to as bilinear.

Thus, equation 10 includes two types of model terms: (a) linear terms (A*K_(i)) and (b) bilinear terms (A*K_(i)).(A*K_(2i)). The star (*) denotes a convolution and the dot (.) denotes a pixel-by-pixel multiplication in the aerial bitmap image.

The bilinear terms (A*K_(i)).(A*K_(2i)) are included to provide a lowest order non-linear response of the resist image to the aerial image. The bilinear terms are the equivalent of a squaring in terms of the frequency support, which results in the doubling of frequencies. Indeed, when two band-limited signals are multiplied, their respective frequency limits are added. In an embodiment, the resist image pixels correspond to half the size of the aerial image pixels. Therefore, it is possible to double the frequency without introducing aliasing errors. With linear terms, the aerial image A is convolved with a kernel, which means that the limit frequency of the aerial image A and the linear term are the same. However, the resist image includes twice the frequency availability. Bilinear convolutions are used to conceptually approximate the reaction-diffusion equations, which are used in the slow and physical models, such as KLA-Tencor's Prolith™ or Brion's LithoCruiser™.

In one implementation, the orthogonal convolution kernels K_(k,n)(r,θ) are selected so that the transformation of the aerial image A into the resist image R has the properties of rotational and mirror symmetry conservation.

The transformation of the aerial image into the resist image has the property of rotational symmetry if the rotation operation does not affect the result of the transformation regardless of whether it is performed before or after the transformation. For example, FIG. 3 shows a transformation that has the property of rotational conservation. As shown in FIG. 3, a resist image 304 that is first calculated with the transformation and then rotated provides the same resist image 306 as a resist image that is calculated by first rotating the aerial image 302 and then applying the transformation to the rotated aerial image 308.

Similarly, the transformation of the aerial image into resist image has the property of mirror symmetry if the symmetry operation does not affect the result of the transformation regardless of whether it is performed before or after the transformation. In other words, if the aerial image A is symmetrical (i.e. A(x,y)=A(−x,y)), then the resist image R is symmetrical too (i.e. R(x,y)=R(−x,y)). For example, FIG. 4 shows a transformation that has the property of symmetry conservation. As shown in FIG. 4, a resist image 404 that is first calculated with the transformation and then subject to a symmetry operation provides the same resist image 408 as a resist image that is calculated by first performing a symmetry operation on the aerial image 402 and then applying the transformation to the mirrored aerial image 404.

In an embodiment, the kernels that provide a rotational conservation for the linear terms are the kernels with n=0 (based on Laguerre polynomials).

In an embodiment, the pair of solutions of convolution kernels that provide a rotational conservation for the bilinear terms are defined as follows: K _(k,n)(r,θ)⁺ =e ^(−r) ² ^(/2) S _(k,n)(r ²)e ^(+inθ) and K _(k′,n)(r,θ)⁻ =e ^(−r) ² ^(/2) S _(k′,n)(r ²)e ^(−inθ), where k,k′, and n are integers.

The radius r can be scaled by the factor r₀ in the linear terms, with r0=√{square root over (2/α)}. As a result, the kernels for the linear terms can be rewritten as follows:

${K_{k,n}\left( {r,\theta} \right)}^{\pm} = {{\mathbb{e}}^{{{- r^{2}}/r_{0}^{2}}/2}{S_{k,n}\left( \frac{r^{2}}{r_{0}^{2}} \right)}{\mathbb{e}}^{{\pm {in}}\;\theta}}$ Similarly, the radius r can be scaled by the factor r₁ in the bilinear terms, with r1=√{square root over (2/α)}. As a result, the kernels for the bilinear terms can be rewritten as follows:

${K_{k,n}\left( {r,\theta} \right)}^{\pm} = {{\mathbb{e}}^{{{- r^{2}}/r_{1}^{2}}/2}{S_{k,n}\left( \frac{r^{2}}{r_{1}^{2}} \right)}{\mathbb{e}}^{{\pm {in}}\;\theta}}$

The proof that it is necessary that the kernels have this form [0065] to follow the property of rotational conservation is evidenced by the following calculation: Indeed, assuming that the two kernels are: K _(k,n)(r,θ)⁺ =e ^(−r) ² ^(/2) S _(k,n)(r ²)e ^(+inθ) and K _(k′,n)(r,θ)⁻ =e ^(−r) ² ^(/2) S _(k′,n′)(r ²)e ^(−in′θ), where k,k′,n and n′ are integers, and assuming that the aerial image A(x,y) is defined as follows A(x, y)=δ(x)δ(y), the bilinear terms become: (A*h ₁).(A*h ₂)=K _(k,n)(r,θ)^(±) .K _(k′,n′)(r,θ)⁼ (A*h ₁).(A*h ₂)=e ^(−r) ² ^(/2) S _(k,n)(r ²)e ^(±inθ) e ^(−r) ² ^(/2) S _(k′,n′)(r ²)e ^(±in′θ) (A*h ₁).(A*h ₂)=e ^(−r) ² S _(k,n)(r ²)S _(k′,n′)(r ²)e ^(±in′θ) e ^(=in′θ) It is necessary that n=−n′.

The proof that it is sufficient that these kernels have this form to follow the property of rotational conservation is evidenced by the following calculation:

For any A: (A*h ₁).(A*h ₂)(x,y)=∫∫A(x−u,y−v)e ^(−r) ² ^(/2) S _(k,n)(r ²)e ^(+inθ) dudv.∫∫A(x−u′,y−v′)e ^(−r′) ² ^(/2) S _(k′,n′)(r′ ²)e ^(−inθ) du′dv′ By first changing the variables from cartesian to polar, the following relations are obtained: x=r₀ cos(θ₀), y=r₀ sin(θ₀), u=r cos(θ), v=r sin(θ), u′=r′ cos(θ′), and v′=r′ sin(θ′). The first integral then becomes: (A*h ₁)=∫∫A(x−u,y−v)e ^(−r) ² ^(/2) S _(k,n)(r ²)e ^(+inθ) dudv If the input image A is rotated by an arbitrary angle α, the new coordinates are given by:

$\begin{bmatrix} {\cos\;\alpha} & {{- \sin}\;\alpha} \\ {\sin\;\alpha} & {\cos\;\alpha} \end{bmatrix}\begin{bmatrix} {x - u} \\ {y - v} \end{bmatrix}$ Then, the first coordinate for the change above is given by:

 = cos  α(x − u) − sin  α(y − v) = cos  α(r₀cos  θ₀ − r cos  θ) − sin  α(r₀sin  θ₀ − r sin  θ) = r₀(cos  θ₀cos  α − sin  θ₀sin  α) − r(cos  θ cos  α − sin  θsinα) = r₀cos (θ₀ + α) − r cos (θ + α) The factor cos(θ₀+α) shows that the output image is rotated. The factor cos(θ+α) by doing a change of variable θ+α→θ, creates a factor e^(−inα) in front of the first integral. The other coordinate of the transform is similar. Similar computations for the second integral are equivalent except that it brings a factor e^(inα), thus cancelling the first factor.

The formalism of section [0058] can be rewritten to ensure that the resist image does not have an imaginary part. For each pair of bilinear terms: z _(n).(A*S _(k) ₁ _(,n)(r ₁ ²)e ^(+inθ) ¹ ).(A*S _(k) ₂ _(,n)(r ₂ ²)e ^(−inθ) ² )+z _(−n).(A*S _(k) ₁ _(,−n)(r ₁ ²)e ^(−inθ) ¹ ).(A*S _(k) ₂ _(,−n)(r ₂ ²)e ^(+inθ) ² ) The above sum can be rewritten as follows: (where

is the real part, ℑ is the imaginary part) d _(n).

└(A*S _(k) ₁ _(,n)(r ₁ ²)e ^(+inθ) ¹ ).(A*S _(k) ₂ _(,n)(r ₂ ²)e ^(−inθ) ² )┘+d _(n)′.ℑ└(A*S _(k) ₁ _(,n)(r ₁ ²)e ^(+inθ) ¹ ).(A*S _(k) ₂ _(,n)(r ₂ ²)e ^(−inθ) ² )┘ Where d _(n) =[z _(n)+ z _(−n) ]/2 and d _(n) ′=[z _(−n)+ z _(n) ]/2 The sum of these two terms is real if and only if both d_(n) and d_(n)′ are real.

The property of symmetry conservation is verified if and only if d_(n)′ is zero. A proof of necessary can be seen by using A(x, y)=δ(x−1)δ(y)+δ(x+1)δ(y) as a counter example. For the proof of sufficiency, the mirror symmetry is computed along the y axis of the resist image: R _(flip)(x,y)=R(−x,y)=(A*h ₁).(A*h ₂)(−x,y) R _(flip)(x,y)=∫∫A(−x−u,y−v)e ^(−r) ² ^(/2) S _(k,n)(r ²)e ^(+inθ) dudv.∫∫A(−x−u′,y−v′)e ^(−r) ² ^(/2) S _(k′,n)(r′ ²)e ^(−inθ′) du′dv′ By a change of variable: u→−u, we have θ→π−θ, du→−du and the integration bounds are switched:

$\begin{matrix} {{R_{flip}\left( {x,y} \right)} =} \\ {{R_{flip}\left( {x,y} \right)} =} \end{matrix}\frac{\begin{matrix} {\int{\int{{A_{flip}\left( {{x - u},{y - v}} \right)}{\mathbb{e}}^{{- r^{2}}/2}{S_{k,n}\left( r^{2} \right)}{\mathbb{e}}^{{- {in}}\;\theta}{\mathbb{d}u}{{\mathbb{d}v} \cdot}}}} \\ {\int{\int{{A_{flip}\left( {{x - u^{\prime}},{y - v^{\prime}}} \right)}{\mathbb{e}}^{{- r^{\prime^{2}}}/2}{S_{k^{\prime},n}\left( r^{\prime 2} \right)}{\mathbb{e}}^{{+ {in}}\;\theta^{\prime}}{\mathbb{d}u^{\prime}}{\mathbb{d}v^{\prime}}}}} \end{matrix}}{\begin{matrix} {\int{\int{{A_{flip}\left( {{x - u},{y - v}} \right)}{\mathbb{e}}^{{- r^{2}}/2}{S_{k,n}\left( r^{2} \right)}{\mathbb{e}}^{{+ {in}}\;\theta}{\mathbb{d}u}{{\mathbb{d}v} \cdot}}}} \\ {\int{\int{{A_{flip}\left( {{x - u^{\prime}},{y - v^{\prime}}} \right)}{\mathbb{e}}^{{- r^{\prime^{2}}}/2}{S_{k^{\prime},n}\left( r^{\prime\; 2} \right)}{\mathbb{e}}^{{- {in}}\;\theta^{\prime}}{\mathbb{d}u^{\prime}}{\mathbb{d}v^{\prime}}}}} \end{matrix}}$ The symmetry condition implies that: R _(flip)(x,y)=∫∫A _(flip)(x−u,y−v)e ^(−r) ² ^(/2) S _(k,n)(r ²)e ^(+inθ) dudv.∫∫A _(flip)(x−u′,y−v′)e ^(−r′) ² ^(/2) S _(k′,n′)(r′ ²)e ^(−inθ′) du′dv′ A number if equal to its complex conjugate if and only if its imaginary part is zero. With the formalism of section [0069], this implies that the coefficient d_(n)′ in front of ℑ└(A*S_(k) ₁ _(,n)(r₁ ²)e^(+inθ) ¹ ).(A*S_(k) ₂ _(,n)(r₂ ²)e^(−inθ) ² )┘ is equal to zero.

Proof for the mirror symmetry along any other axes is done by combining the property of rotation conservation of section [0068] and the property of mirror symetry conservation along the y-axis.

An embodiment of the invention provides a fast implementation of the model after it has been fitted. In equation 10, the two kernels K_(i) and K_(2i) of the bilinear terms can be different. An eigenvalue transformation of the kernels allows one to implement only transformed terms where K_(i) and K_(2i) are the same. This reduces the number of kernels to be applied, from n² pairs to n terms (on average). In the formalism that follows, the convolutions are treated as 1D for simplicity of the writing, but by simple bookkeeping, the procedure is extended to 2D filters. The procedure is as follows:

For a given, fixed n, the following equation is computed:

$R = {\sum\limits_{k,k^{\prime}}{{z_{k,k^{\prime}}\left( {A*K_{k,n}} \right)} \cdot \left( {A*K_{k^{\prime},n}} \right)}}$ The matrix Z=└z_(k,k′)┘ is built and its eigenvalues w_(i) and corresponding eigenvectors V_(i) are computed. The matrix K=└K_(k,n)┘ is created and the vectors Y_(i)=K. V _(i) are computed. The following relation is obtained

$R = {\sum\limits_{i}{w_{i}{{{A*Y_{i}}}^{2}.}}}$ The number of convolutions needed increases linearly, instead of square. The procedure is repeated for every n. This computation may be done once and applied to various aerial images.

Referring now to FIG. 5, this figure is an exemplary flowchart illustrating a process 500 of obtaining a resist image R in accordance with an embodiment of the invention. The process 500 begins at procedure 505 where an aerial image A of the mask pattern is determined at substrate level. Calculation of the aerial image A may be performed in the same manner as in procedure 205 of process 200. After obtaining the aerial image A, the process 500 proceeds to procedure 510 where the plurality of fitting coefficients c_(i) and d_(i) of equation 10 are selected. Fitting coefficients c_(i) and d_(i) are linear coefficients that are iteratively adjusted to determine the resist image. Fitting coefficients c_(i),d_(i), r₀ and r₁ can be automatically generated or selected by the user. Subsequently, a resist image R is calculated using equation 10 (procedure 515). The result of the transformation of the aerial image A by equation 10 provides a resist image R that is converted into a critical dimension (CD) value at procedure 520. The CD value (or contour) obtained is compared with a target value at procedure 520. The target value is obtained by measuring an actual substrate. The actual substrate to be measured is prepared by use of the optical tools which were used to obtain the aerial image. Then, the method proceeds to procedure 530 where a determination is made as to whether the CD value matches the target CD value. When the CD value is within a predetermined range, i.e., equal to or close to the target value (step 18: yes), the process ends (procedure 535) and the CD value represents an accurate prediction of the critical dimension of the mask pattern in the resist. On the other hand, when the CD value is not within the predetermined range, the method proceeds back to procedure 510 where new values of fitting coefficients c_(i), d_(i), r₀ and r₁ are selected.

The iterative process of procedures 510-530 can be performed using different types of regression schemes. The fitting of coefficients c_(i) and d_(i) consists of solving the over-determined system Bx=C, where x represents a vector that contains all the coefficients c_(i) and d_(i) to be determined, B represents a matrix where each column corresponds to a resist term and each row corresponds to one gauge (i.e. the term B_(i,j) is the intensity of term image j at jauge i) and the vector C contains the target intensity. Since the system is over-determined, it cannot be solved exactly. Therefore, in one implementation, the system is solved by determining the coefficients of the vector x that minimizes ∥Ax−B∥²+α²∥x∥² (Eq. 11). The term α is fixed and depends on the amount of noise. In the context of model fitting, this approach is known as a ridge regression, a Tikhonov regularization or Philips regularization. The parameter α can be obtained by using a Wahba routine, as will be appreciated by one skilled in the art. Once the parameter α is determined, equation 11 can be solved to obtain the fitting coefficients.

The use of the process shown in the embodiments of FIGS. 2 and 5 for calculating the resist image provides significant benefits over conventional approaches. Contrary to conventional approaches, the kernels used in embodiments of the invention form a complete basis. This allows one to improve the fitting accuracy, as it is possible to cover the space of all possible images.

Moreover, contrary to conventional approaches, the convolution kernels form an orthogonal basis. The orthogonality lowers the risk cancellations of terms in Equation 10 during the fitting process, due to degeneracy or redundancy. As a result, it is possible to cover more of the vector space with the same number of kernels.

In addition, it will be appreciated that the process of FIGS. 2 and 5 can be carried out with a much more reduced set of parameters and thus requires less user interaction as compared to a conventional simulation method. For example, the user does not have to experiment several search ranges nor does the user have to turn on/off several terms. Instead, the search is straightforward and there are a limited number of settings.

Further, it will be appreciated that the resist image R can be determined with a limited number of terms in Equation 10. For example, the use of the following collection of vector basis yields good results:

TABLE 1 Type n k1 K2 Radius of Gaussian Linear (HL) 0 0 r₀ Linear (HL) 0 1 r₀ Linear (HL) 0 2 r₀ Linear (HL) 0 3 r₀ Linear (HL) 0 4 r₀ Bilinear (HB) 0 0 0 r₁ Bilinear (HB) 0 0 1 r₁ Bilinear (HB) 0 0 2 r₁ Bilinear (HB) 0 1 1 r₁ Bilinear (HB) 0 1 2 r₁ Bilinear (HB) 0 2 2 r₁ Bilinear (HB) 1 1 1 r₁ Bilinear (HB) 1 1 2 r₁ Bilinear (HB) 1 2 2 r₁ where n, k1 and k2 are the coefficients of the solution kernels

${K_{k,n}\left( {r,\theta} \right)}^{\pm} = {{\mathbb{e}}^{- \frac{r^{2}}{r_{1}^{2}/2}}{S_{k,n}\left( \frac{r^{2}}{r_{1}^{2}} \right)}{\mathbb{e}}^{{\pm {in}}\;\theta}}$ for a bilinear term, and

${K_{k,n}\left( {r,\theta} \right)}^{\pm} = {{\mathbb{e}}^{- \frac{r^{2}}{r_{0}^{2}/2}}{S_{k,0}\left( \frac{r^{2}}{r_{0}^{2}} \right)}}$ for a linear term.

In an embodiment, the search range of radii is defined as follows:

minimum maximum step r₀ 20 100 10 r₁ 20 100 10

The calculated resist image R can be used to determine the critical dimension of a pattern, the contour of a pattern and whether some contours of resist patterns are merged or cut. The calculated image can also be used to determine the size of optical proximity correction features. Further, it will be appreciated that the procedures and methods defined hereinabove can also be used to calculate an etch image.

It will be appreciated that the different acts involved in calculating the image of the mask in the resist may be executed according to machine executable instructions or codes. These machine executable instructions may be embedded in a data storage medium. A processor may be configured to execute the instructions. In one embodiment, the resist image computations are carried out in FPGA (Field Programmable Gate Array) in a similar manner as the aerial image computations.

Software functionalities of a computer system involving programming, including executable codes, may be used to implement the above described imaging model. The software code may be executable by a general-purpose computer. In operation, the code and possibly the associated data records may be stored within a general-purpose computer platform. At other times, however, the software may be stored at other locations and/or transported for loading into an appropriate general-purpose computer system. Hence, the embodiments discussed above involve one or more software or computer products in the form of one or more modules of code carried by at least one machine-readable medium. Execution of such codes by a processor of the computer system enables the platform to implement the functions in essentially the manner performed in the embodiments discussed and illustrated herein.

As used herein, terms such as computer or machine “readable medium” refer to any medium that participates in providing instructions to a processor for execution. Such a medium may take many forms, including but not limited to, non-volatile media, volatile media, and transmission media. Non-volatile media include, for example, optical or magnetic disks, such as any of the storage devices in any computer(s) operating as discussed above. Volatile media include dynamic memory, such as the main memory of a computer system. Physical transmission media include coaxial cables, copper wires and fiber optics, including the wires that comprise a bus within a computer system. Carrier-wave transmission media can take the form of electric or electromagnetic signals, or acoustic or light waves such as those generated during radio frequency (RF) and infrared (IR) data communications. Common forms of computer-readable media therefore include, for example: a floppy disk, a flexible disk, hard disk, magnetic tape, any other magnetic medium, a CD-ROM, DVD, any other optical medium, less commonly used media such as punch cards, paper tape, any other physical medium with patterns of holes, a RAM, a PROM, and EPROM, a FLASH-EPROM, any other memory chip or cartridge, a carrier wave transporting data or instructions, cables or links transporting such a carrier wave, or any other medium from which a computer can read or send programming codes and/or data. Many of these forms of computer readable media may be involved in carrying one or more sequences of one or more instructions to a processor for execution.

While specific embodiments of the invention have been described above, it will be appreciated that the invention may be practiced other than as described. The description is not intended to limit the invention.

Although specific reference may be made in this text to the use of lithographic apparatus in the manufacture of ICs, it should be understood that the lithographic apparatus described herein may have other applications, such as the manufacture of integrated optical systems, guidance and detection patterns for magnetic domain memories, flat-panel displays, liquid-crystal displays (LCDs), thin-film magnetic heads, etc. The skilled artisan will appreciate that, in the context of such alternative applications, any use of the terms “wafer” or “die” herein may be considered as synonymous with the more general terms “substrate” or “target portion,” respectively. The substrate referred to herein may be processed, before or after exposure, in, for example, a track (a tool that typically applies a layer of resist to a substrate and develops the exposed resist), a metrology tool and/or an inspection tool. Where applicable, the disclosure herein may be applied to such and other substrate processing tools. Further, the substrate may be processed more than once, for example in order to create a multi-layer IC, so that the term substrate used herein may also refer to a substrate that already contains multiple processed layers.

The lithographic apparatus may also be of a type wherein a surface of the substrate is immersed in a liquid having a relatively high refractive index, e.g., water, so as to fill a space between a final element of the projection system and the substrate Immersion liquids may also be applied to other spaces in the lithographic apparatus, for example between the patterning device and a first element of the projection system Immersion techniques are well known in the art for increasing the numerical aperture of projection systems.

The terms “radiation” and “beam” used herein encompass all types of electromagnetic radiation, including ultraviolet (UV) radiation (e.g., having a wavelength of or about 365, 248, 193, 157 or 126 nm) and extreme ultra-violet (EUV) radiation (e.g., having a wavelength in the range of 5-20 nm), as well as particle beams, such as ion beams or electron beams.

The term “lens,” where the context allows, may refer to any one or combination of various types of optical components, including refractive, reflective, magnetic, electromagnetic and electrostatic optical components.

While specific embodiments of the invention have been described above, it will be appreciated that the invention may be practiced other ways than as described. For example, the invention may take the form of a computer program containing one or more sequences of machine-readable instructions describing a method as disclosed above, or a data storage medium (e.g., semiconductor memory, magnetic or optical disk) having such a computer program stored therein.

The descriptions above are intended to be illustrative, not limiting. Thus, it will be apparent to those skilled in the art that modifications may be made to the invention as described without departing from the scope of the claims set out below. 

What is claimed is:
 1. A method implemented by a computer processor for simulating a resist image, wherein the resist image comprises an image of a portion of a design layout to be lithographically projected onto a resist coated on a substrate, the method comprising: determining an aerial image of the portion of the design layout at substrate level; and transforming the aerial image into the resist image by convolving the aerial image with selected convolution kernels having both linear and bilinear terms.
 2. The method of claim 1, wherein the transforming step includes: selecting the convolution kernels that form a complete basis, covering a vector space of all possible resist images.
 3. The method of claim 1, wherein the transforming step includes: selecting the convolution kernels that form an orthogonal basis, wherein correlation between the selected convolution kernels is minimal.
 4. The method of claim 3, wherein the selected convolution kernels comprise solutions of a two-dimensional quantum harmonic oscillator.
 5. The method of claim 4, wherein the solutions are based on Laguerre and Whittaker functions that modulate Gaussian kernels and form the orthogonal basis.
 6. The method of claim 1, wherein the transforming step includes: selecting the convolution kernels to preserve mirror symmetry of the resist image.
 7. The method of claim 1, wherein the transforming step further includes: selecting the convolution kernels such that a pair of convolution kernels has rotationally-independent bilinear terms, thereby preserving rotational symmetry of the resist image.
 8. The method of claim 1, wherein the simulation of the resist image is based on an empirical model.
 9. The method of claim 8, wherein the bilinear terms in the convolution kernels conceptually approximate chemical reaction-diffusion equations in the empirical model.
 10. A non-transitory computer program product having computer executable instructions for a computer processor to perform a method for simulating a resist image, wherein the resist image comprises an image of a portion of a design layout to be lithographically projected onto a resist coated on a substrate, the method comprising: determining an aerial image of the portion of the design layout at substrate level; and transforming the aerial image into the resist image by convolving the aerial image with selected convolution kernels having both linear and bilinear terms.
 11. The computer program product of claim 10, wherein the transforming step includes: selecting the convolution kernels that form a complete basis, covering a vector space of all possible resist images.
 12. The computer program product of claim 10, wherein the transforming step includes: selecting the convolution kernels that form an orthogonal basis, wherein correlation between the selected convolution kernels is minimal.
 13. The computer program product of claim 12, wherein the selected convolution kernels comprise solutions of a two-dimensional quantum harmonic oscillator.
 14. The computer program product of claim 13, wherein the solutions are based on Laguerre and Whittaker functions that modulate Gaussian kernels and form the orthogonal basis.
 15. The computer program product of claim 10, wherein the transforming step includes: selecting the convolution kernels to preserve mirror and rotational symmetry of the resist image.
 16. The computer program product of claim 15, wherein the transforming step further includes: selecting the convolution kernels such that a pair of convolution kernels has rotationally-independent bilinear terms, thereby preserving rotational symmetry of the resist image.
 17. The computer program product of claim 10, wherein the simulation of the resist image is based on an empirical model.
 18. The computer program product of claim 17, wherein the bilinear terms in the convolution kernels conceptually approximate chemical reaction-diffusion equations in the empirical model.
 19. A method implemented by a computer processor for simulating an image of a portion of a design layout to be lithographically etched on a substrate, the method comprising: determining an aerial image of the portion of the design layout at substrate level; and transforming the aerial image into an etch image by convolving the aerial image with selected convolution kernels, having both linear and bilinear terms, wherein the etch image is representative of the image of the design layout lithographically etched on the substrate.
 20. The method of claim 19, wherein the transforming step includes: selecting the convolution kernels that form a complete basis, covering a vector space of all possible etch images.
 21. The method of claim 19, wherein the transforming step includes: selecting the convolution kernels that form an orthogonal basis, wherein correlation between the selected convolution kernels is minimal.
 22. The method of claim 21, wherein the selected convolution kernels comprise solutions of a two-dimensional quantum harmonic oscillator.
 23. The method of claim 22, wherein the solutions are based on Laguerre and Whittaker functions that modulate Gaussian kernels and form the orthogonal basis.
 24. The method of claim 19, wherein the transforming step includes: selecting the convolution kernels to preserve mirror symmetry of the etch image.
 25. The method of claim 19, wherein the transforming step further includes: selecting the convolution kernels such that a pair of convolution kernels has rotationally-independent bilinear terms, thereby preserving rotational symmetry of the etch image.
 26. The method of claim 19, wherein the simulation of the etch image is based on an empirical model. 